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Abstract 

In this paper a library for spin-angular integration in LS-coupling for many-electron atoms is pre- 
sented. The software is an implementation of a methodology based on the second quantization in coupled 
tensorial form, the angular momentum theory in 3 spaces (orbital, spin and quasispin), and the graphical 
technique of angular momentum. This implementation extends applications of the relevant methodology 
to open /- shells and leads to faster execution of angular integration codes. The possibility of using some 
library routines for solving various angular momentum problems in atomic physics is also discussed. 
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PROGRAM SUMMARY 



Title of program : SAI Catalogue identifier: 

Program obtainable from: Institute of Theoretical Physics and Astronomy, A. Gostauto 12, 2600 Vilnius, LITHUA- 
NIA. 

E-mail: gaigalas@itpa.lt 

Computer for which the library is designed and others on which it has been tested: 
Computers: Pentium-based PCs; 

Installations: Institute of Theoretical Physics and Astronomy, A. Gostauto 12, 2600 Vilnius, LITHUANIA. 
Programming language used in the new version: FORTRAN 77 [1] 
Peripherals used: terminal, disk 
No. of bits in a word: 32 
No. of processors used: 1 

Has the code been vectorised or parallelized?: no 

No. of bytes in distributed program, including test data, etc. : 442 297 bytes 
Distribution format: compressed tar file 

Additional keywords : atomic structure, configuration interaction, reduced coefficients of fractional parentage, irre- 
ducible tensors, angular momentum theory in three spaces (orbital, spin and quasispin), second quantization in the 
coupled tcnsorial form, recoupling coefficients, Slater integrals, complex atom, correlation, wave functions, bound 
states, I/S-coupling, /- shell states. 

Nature of physical problem 

Theoretical determination of atomic orbitals, energy levels and radiative transition data requires the calculation of 
matrix elements of the relevant physical operators (see the multiconfiguration Hartree-Fock method [2], for example). 
The matrix elements of arbitrary operator can generally be expressed as 
coef(i,j) {^USiWT^ULiSj), 

where T (fc;fcs) is a tensorial operator of ranks ki, k 3 . The program calculates the spin-angular part for matrix 
elements {'y i LiSi\\T < - klks) \\'yjLjSj) of one- and/or two-particle operator T (felfcs) . 

Method of solution 

This program is created involving the angular methodology of [3-6] . It has been extended to include partially filled 
/- shells in wave function expansions. The classification of terms is identical to that described in [6]. 

Restrictions on the complexity of the problem 
The non-orthogonal orbitals are not supported. 

References 

[1] Microsoft, Fortran Powerstation Programmer's Guide (Microsoft Corporation, 1995). 

[2] C. Froese Fischer, T. Brage and P. Jonsson, Computational Atomic Structure. An MCHF Approach (Institute of 

Physics Publishing, Bristol/Philadelphia, 1997). 
[3] G. Gaigalas and Z. Rudzikas, J. Phys. B: At. Mol. Phys. 29 (1996) 3303. 
[4] G. Gaigalas, Z. Rudzikas and C. Froese Fischer, J. Phys. B: At. Mol. Phys. 30 (1997) 3747. 
[5] G. Gaigalas, A. Bernotas, Z. Rudzikas and C. Froese Fischer, Physica Scripta 57 (1998) 207. 
[6] G. Gaigalas, Z. Rudzikas and C. Froese Fischer, Atomic Data and Nuclear Data Tables 70 (1998) 1. 



2 



LONG WRITE-UP 



1 Introduction 

Models of many-electron atoms and ions require both relativistic and correlation effects to be taken into 
account for getting very precise characteristics of atoms and ions. This can be done, for example, by using 
multiconfiguration Haretree-Fock, configuration interaction method, various versions of perturbation theory, 
or semiempirical methods [1] . All of them require the calculation of matrix elements of physical operators or 
effective operators from perturbation theory. The symmetry properties of atomic states allow the calculation 
of matrix elements to be divided into the calculation of spin-angular terms and the accompanying radial 
integrals. The latter are the more straightforward and can be handled by methods such as used in the MCHF 
atomic structure package (ATSP_MCHF) [2, 3]. Such packages have a modular structure, and the modules for 
calculation of the spin-angular part for matrix elements of any operator can easily be replaced with new 
one, which is more efficient for large scale computation of open shell atoms. 

This paper describes modules of such a sort. The new program (library for integration over spin-angular 
variables in atomic theory - library SAI) is based on combination of the angular momentum theory, on 
the concept of irreducible tensorial sets, on a generalized graphical approach, on the second quantization in 
coupled tensorial form, on the quasispin approach and on standard quantities like the reduced coefficients 
of fractional parentage, completely reduced matrix elements of the unit tensors as well as a number of other 
completely reduced matrix elements occurring in various products of electron creation and annihilation 
operators [4]. The program for calculation of standard quantities have been published as library SQ [5]. The 
program presented in this paper uses the library SQ for evaluation of these standard quantities. The new 
module for the calculation of spin-angular coefficients is faster compared to previous ones and can treat 
configurations with open f N shell. 

The theoretical background is presented in Section 2, program organization is outlined in Section 3, the 
description of the use of this library for other programs is given in Section 4. The program installation is 
presented in Section 5. Finally, a few examples are given in Section 6. 

2 Notations and Methodology of Angular Integrations 

2.1 Matrix Elements Between Complex Configurations 

According to the approach [4, 6], a general expression of a submatrix element for any two-particle operator 
between functions with u open shells, can be written as follows: 

(V„ {LS)\\GU U (L'S')) 

mli.ndj .n'.l'n'.V. 

E EE (-ife'M^A^A^A^s) 

xT (nA, njXjXK, ^;,A fc ™, A ket , 5, r) R (X u A,, Aj, X' v A 6 ™, A ket , r) , (1) 

where A = l,s, A\ ra = (L u L h L\, L' j ) bra , A b s ra = (S u Sj, S' { , srf™ and T refers to the array of cou- 
pling parameters connecting the recoupling matrix R (Aj, Xj, A-, X'j,A bra , A ket , T) to the submatrix element 
T ( y riiXi 7 njXj 7 n' i X' i7 n' : jXj,A bra 7 A ket } E,r). The expression (1) has summation over intermediate ranks K\2, 
a 12 , k' 12 , a' 12 , K u K s in T (n i X i ,n : jX : j,n' i X' i ,n' j X' : j,A hra ,A ket ,E,T). 

So, to calculate the spin-angular part of a submatrix element of this type, one has to obtain: 

1. Recoupling matrix R (Aj, Aj, A^, X' 3 A bra 1 A ket 1 T) . This recoupling matrix accounts for the change in 
going from matrix element (ip u (LS) \ \G\ \ip u (L'S')), which has u open shells in the bra and ket functions, 
to the submatrix element T (rijAj, rij Xj ;, n-A^, n'j X'j, A bra , A fce *, S, r) , which has only the shells being 
acted upon by the two-particle operator in its bra and ket functions. 
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2. Submatrix elements T (rijAj, njXj,n' i X' i ,n'jX' j ,A bra , A ket ,E, T). 

3. Phase factor A. 

4. 6' (njAj,njAj,n-A^,rijA^,5), which is proportional to the radial part. It consists of a submatrix element 
(ni\inj\j\\g( KlK2k, ' 7l ' T2k ^\\ip u (L'S')), and in some cases of simple factors and 3nj-coefficients (for more 
details see [6]). 

Some important points to note are the following: 

1 . The recoupling matrices R(X i ,X j ,X' i , X'j , A bra ,A ket , T) in our approach are much simpler than in other 
known approaches. We have obtained their analytical expressions in terms of just 6j- and 9j-cocfficients. 
That is why we choose a special form of operator in second quantization, where second quantization operators 
acting upon the same shell are tensorially coupled together. 

2. The tensorial part of a two-particle operator is expressed in terms of (products of) operators of the 
type A( kk t (nA, ~), (nA, 5), C^ kk \nX, 5), D^ ls \ E^ kk \nX,E). Their explicit expressions are (2)-(6): 



,(</A) 
rn ql 



,(?A) 

"1,2 



(KlCl) 



<4f ? x 



(kicti) 



<x 5 >< <l 



(KlCTl) 



x a. 



(</A) 
™ 9 3 



(K2CT2) 



(«2C 2 ) 



m ql m q2 



(KlCTl) 



a^ A , } x a<f ] 



(K2CT2) 



(fcfc) 



(2) 
(3) 

(4) 
(5) 
(6) 



We denote their submatrix elements by T (rijAj, n 3 Xj, n-A-, ; Aj •, A bra , A ket , 5, r) . The parameter T repre- 
sents the whole array of parameters connecting the recoupling matrix R (Aj, Aj, A-, X'j,A bra , A ket , r) to the 
submatrix element T (n^A^, n^Aj, n^AJ, n^-A^ , A hra , A fce *, S, r) . It is worth noting that each of the tensorial 
quantities (2)-(6) act upon one and the same shell. So, all the advantages of tensor algebra and the quasispin 
formalism may be efficiently exploited in the process of their calculation. 

We obtain the submatrix elements of operator (2) by using straightforwardly the Wigner-Eckart theorem 
in quasispin space: 



(l N aQLS\\a^\\l N ' a'Q'L'S') 



[Q] 



-1/2 



Q' 1/2 Q 
M'q m q M Q 



(l aQLS\\\a 



(qls) I 



a'Q'L'S 1 



(7) 



where the last multiplier in (7) is the so-called reduced coefficient of fractional parentage (RCFP) and we 
use a shorthand notation (2k + 1) • ... = [k, ...]. The details of the use of quasispin approach are discussed in 
monograph [1]. 

The value of a submatrix element of operator (3) is obtained by basing ourselves on (33), (34) in [7]. In 
the other three cases (4), (5), (6) we obtain them by using (5.16) of Rudzikas [1]: 

r -1 (kk) 

{nl N aQLS\\ F (Klffl) (nA) x G (K2CTa) (nA) \\nl N ' a'Q'L'S') 



= (-1) 



L+S+L' + S'+2k 



[k] {nl N aQLS\\F^°^ [nX)\\nl N " a"Q"L"S") 



a"Q"L"S" 



x(nl N " a" ■Q"L"S"\\G^(n\) \\nl N ' a'Q'L'S') ^ *„ } { J 5 5"}' 



(8) 
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where F { ^ a ^ (nX) and G {k ^ (nX) is one of (2) or (3) and the submatrix elements correspondingly are 
defined by (7) and (33), (34) in [7]. N" is defined by second quantization operators occurring in 
and G< K2<T2 ) (nX). 

As is seen, by using this approach, the calculation of the angular parts of matrix elements between 
functions with u open shells ends up as a calculation of submatrix elements of tensors (2), (3) within single 
shell of equivalent electrons. As these completely reduced submatrix elements (reduced in the quasispin, 
orbital and spin spaces) do not depend on the occupation number of the shell, the tables for them arc 
reduced considerably in comparison with the tables of analogous submatrix elements of tensorial quantities 
U , V klk2 (Jucys and Savukynas [8] or Cowan [9]) and the tables of fractional parentage coefficients (CFP). 
That is why the expressions obtained are very useful in practical calculations. This is extremely important 
for the /- shell, where the number of CFPs for f 1 - / 14 equals 54408 [10], whereas the number of RCFPs, 
taking into account the transposition symmetry property of RCFP, is only 14161 - of which only 3624 are 
nonzero [11]. 

We do not present details on obtaining phase factors A and 6' (rijAj, rijXj, n-A^, n'jX'j, S) , since no essential 
generalizations may be made here; those are possible only after a particular operator is chosen (for more 
details see [6, 12]). 

2.2 The Electrostatic Electron Interaction, Spin-Spin and Spin-Other-Orbit 
Operators 

The electrostatic (Coulomb) electron interaction operator H Coulomb itself contains the tensorial structure 



rrCoulomb _ rr(/c/c0,000) 
n — Coulomb 



(9) 



and its submatrix clement is: 



UiXiUjXj 



H. 



(fcfcO,000) 
Coulomb 



n v \i.nj>\j^ = 2[fc] 1/2 [l, C [k) (l 3 C [k) i? fe (n^rv^ ,n 3 l 3 n 3 d r ) . (10) 



The spin-spin operator H ss itself contains tensorial structure of two different types, summed over k: 



H ss = J2[ H ss + lh ~ 12,112) +H ss~ 



lfe+12,112) 



(11) 



Their submatrix elements are: 



(riiXiUj^ 



^j(fe+lfe-12,112) 



niiXiinyXji^ 



= -^(2fc + 3) (5) (*i||c (fc+1) ||*<') [ijjc^jl^ N k - X (n^n^^nek^d,,) , (12) 



{niXiUjXj 



jy(fe-lfc+12,112) 



Xi'riji Xji^j 



A^/(2fc + 3 )< 6 > {h 



c (k-i) 



) (Zj C (fe+1) l r ^N k - 1 (njljmlunjdj.nidi,) , (13) 



where we use a shorthand notation (2k + 3)^ = (2k + 3) (2k + 2) (2k + 1) (2k) (2k — 1) and radial integral 
(12), (13) is defined as in Glass and Hibbcrt [13]: 



2 roo roc fc 

) = X / Pi (n) Pj (r 2 ) e(n - r 2 )P t , (n) Py (r 2 ) dr x dr 2 

4 Jo Jo r 1 



(14) 



where e(x) is a Heaviside step-function, 



e(x) = 



1; for x > 0, 
0; for x < 0. 



(15) 
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The spin-other-orbit operator H so ° itself contains tensorial structure of six different types, summed over 



TTSso — [rWfc-lfcl.lOl) , rj-(fc-lfcl,011) , rx(fcfcl,101) , rj-(fefel,011) , rx(/c+l/cl,101) , rj-(fe+lfel,011) 

n / j ^sso ' ^sso ' ^sso ' ^sso ' ^sso ' ^sso 



(16) 



Their submatrix elements are: 



(jiiXirijXj 



jj(k — lkl,<Tl(T2l) 



niiXiirijtXjt^j 



= 2 • 2 a2 {(2k - 1) (2k + 1) (h + h> - k + 1) (fc - h + U,) (k + h- l v ) (k + h + l v + 1)} 1/2 
x(ky 1,2 (l l C (fe) k^{lj C (fe) Ij^j N k ~ 2 (njlimlunylymdv) , (17) 



(niXtHjXj 



Jj(k+lkl,(7i(T2l) 



n^X^rijiXj^j 



= 2 • 2 CT2 {(2fc + 1) (2fc + 3) (h + l v -k)(k-h + l v + l)(k + h- l v +l)(k + k + l v + 2)} 1/2 
x(fc + l)~ 1/2 (^ C {k) U')(lj C (fe) Nk {n i l i njlj,n i ,l i ,nj l l j ,). 



(18) 



{niXiUjXj 



H 



{kkl,a\ (72 1) 



T^i' ^i' Tlj' 



= -2 ■ 2 a2 (2k + 1) 1/2 (l, \c {k) \ k>) (lj \c {k) II l r ^j {(fc (fc + I)) -172 

x (h (h + 1) - fc(fc + 1) - l v (l v + 1)) {(fc + l)iV fe - 2 (njljnili,n r l r nidi,) 

— kN k (nilinjlj,riiiliinjilji)} — 2 (fc (fc + l)) 1 ^ 2 V k ~ x (nilinjlj,riiiliinjilji)^ . 

The radial integrals in (17)-(19) are defined by (14) and below (see Glass and Hibbert [13]): 

a 2 f°° f°° r k ~ x d 

V (nilinjlj^i'li'rij'lf) = — / Pi (n) Pj (r-2) -r+2 r ^7T~ Pi ' ( ri ) p i' ( r2 ) dndr 2 . 

4 Jo Jo r > or\ 



(19) 



(20) 



Now we have all we need (the operators with tensorial structure and their submatrix elements) for 
obtaining the values of a matrix element of these operators for any number of open shells in bra and ket 
functions. This lets us exploit all advantages of the approach by [6]. 

The spin-spin and spin-other-orbit operators themselves generally contain tensorial structure of several 
different types. Therefore the expression (1) must be used separately for each possible tensorial structure 
for performing spin-angular integrations according to [6] . Each type of tensorial structure is associated with 
a different type of recoupling matrix R (X il Xj, X' t , X' 3 ., A 6 ™, A ket , r) and with different matrix elements of 
standard tensorial quantities T (riiXi, rijXj, n\X\, n'^X'^ ■, A bra , A ket , S, r) . 

The one-particle operators are treated in a similar manner. Their expressions are much simpler and 
therefore we do not present them here, for brevity They may be found in [14]. 



3 Description of the Library 

The library SAI presented in this paper is aimed at the spin-angular integration for any one- and two- 
particle operator. It is a separate unit and can easily be adapted to existing codes such as the MCHF atomic 
structure package (ATSP_MCHF) [2, 3] or can easily be used to create a new one. It contains four modules - 
SAI_RECLS, SAI_SQ, SAIJJORE and SAI_DUDU. They are classified according to the methodology presented in 
papers [4, 6, 7, 12], and adhere to the principles of modular programming (although FORTRAN 77 does not 
fully support this). 
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The module SAI_SQ was published separately [5]. In this paper we will discuss briefly the new modules 
and subroutines contained therein. The author does not attempt to describe in detail all the subroutines 
belonging to separate libraries. In order to give the reader a more complete view of the implementation 
of methodology published in papers [4, 6, 7, 11, 12, 15], only the main subroutines are described. Special 
attention is given to the description of subroutines that other programs must call for use of this library. The 
subroutines described in more detail can easily be used separately from the complete library SAI, too. These 
subroutines are useful for creating new programs or modifying old ones, even those based on the traditional 
methodology of angular calculations described by Fano [16]. 

3.1 SAI_RECLS 

This library contains 20 routines for calculation of recoupling matrices 

R A,, X'^X'j, A 6 ™, A fcet , T) = R {k, I'ntiT, Af et , r,) R (a, s, s, s, A b s ra , A ket ,T s ) . (21) 

For more details see [6] (Section 4) . Most of the subroutines from this module use common blocks CONSTS 
and MEDEFN from ATSP_MCHF [2, 3]. 

DLSAn, n=l, 2, ... 6, evaluates respectively coefficients C\, C5, C 2 , C4, C3, C 6 , defined in equations 
(15), (16), (23), (21), (17), (25) of [6]. 

RECOUPO checks the angular momentum selection rules for the recoupling coefficients. For example, 
it uses the expression (18) in one interacting shell case (see [6]). 



RECOUPn, n=2, 3, 4, evaluates the Kronecker delta functions or calculates the recoupling coefficients 
respectively for the scalar operator 



A™ ( ni h) x (n 2 l 2 ) 



(o) 



(22) 



A^ (mh) x (rijlj) x (n m l m ) 



(fe) 



(0) 



(23) 



(*0 



A^ (mh) x B^ (n 2 h) x C (fe3) {n 3 l 3 ) x {n 4 l 4 ) 



(*0 



(0) 



(24) 



defined in equations (22), (26), (33) of [6]. The A^ k \ B {k2 \ and may be simple or composite 

tensorial operators. 

RLSPO evaluates the Kronecker delta functions 5(Li,L'A for one and two interacting shells (see (14) 
and (19) in [6]). 

RLSPOO evaluates the Kronecker delta functions 5 (Li, L-) for three and four interacting shells (see (24) 
and (27) in [6]). 



RLSPn, n=l, 2, 3 evaluates the Kronecker delta functions or calculates the recoupling coefficients 
respectively for the non-scalar operator 



A {k) (mh).. 



(25) 



A {k ^ (mil) x B {k ^ (n 2 l 2 ) , 



(fe) 



(26) 



A^ (mh) x B^ (ml,) x C {ki) (n m l m ) 



(fc 3 ) 



(fe) 



(27) 
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defined in equations (14), (19), (24) of [6]. This routine evaluates the Kronecker delta functions or calculates 
respectively the first or second part of recoupling coefficients for the non-scalar operator 

[A^ (mil) x {n 2 l 2 ) 
in the case n=4a, 4b. 



C (fc4) {n 3 l 3 ) x £>(**> (n 4 l 4 



(28) 



3.2 SAI_NORE 



This library is for calculating the spin angular parts of matrix elements for a scalar two-particle operator. It 
contains 18 subroutines. Most of the subroutines from this module use common blocks CONSTS and MEDEFN 
from ATSP_MCHF [2, 3]. 

From (1) we see that the matrix element of any two-particle operator can be written as a sum over 
all possible sets of active shell quantum numbers rnk. The systematic analysis of [6] aims to minimize the 
number of distributions, which is necessary to obtain the matrix elements of any two-electron operator, 
when the bra and ket functions consist of arbitrary number of shells. Table 1 lists all these distributions and 
at the same time the expressions used by each of the subroutines NONRELAT1, NONRELAT2, NON- 
RELAT31, NONRELAT32, NONRELAT33, NONREL AT41 , NONRELAT51, NONRELAT52, 
NONRELAT53. The numbering of expressions is the same as in paper [6], where all these expressions 
are presented. As the structure of all the subroutines mentioned earlier is the same, and only different 
expressions are used and different subroutines are called, we will discuss in more detail only one of these 
subroutines. 

NONRELAT1 is meant for finding spin-angular coefficients for the distributions aaaa, af3a(3 and 
(3 a Pa. 

In the aaaa case, the program uses expression (5) from [6]. In this case 

G (-0 ~ E E ( nA ' S ) A « 0) ( nA ' S ) + e "«- ( nA ' S ) A (SI-p ( " A ' S) ] ' (29) 

Kl2,0'12,K' 1 2!0-! 2 P 



where 



4° 0) («A,S) 



a<'">xa<'-> ^x^xa^) 10 ^ 



(00) 



(30) 



and 



4S,->A,S) = 



a (l a s) x ~(Z qS ) 



(fefe) 
P,-P 



(31) 



and 



©//a (n\, S) = 9// a (n a X a , n a X a , n a \ a , n a \ a , S) 

= l -{-lf- p [K 1 ,a U K 2 ,a 2 }- 1/2 (n a X a n a X a \\g^ K2k ^V\\n a X a n a X a ) 

©776 (nX, S) = 9//6 (n a X a , n a X a , n a X a , n a X a , S) 

= (-l) fe - p+1 (n a A a n a A a || 5 (— fe '— fc )||n a A a n a A a ){ * ^ M { ^ ^ J 



(32) 



(33) 



The value of the reduced matrix element of operator (30) is found by subroutine WWLS1, and that of 
(31) by Wl (see module SAI_SQ in [5]). The value of coefficient Q IIa (nX, S) is calculated by COULOMBLS 
subroutine, because this coefficient, to the accuracy of a factor and a phase, is equal to that part of the two- 
electron submatrix elements of electrostatic interaction, which this subroutine is calculating. The coefficient 
@776 (nX, S) is found by COULOMBLS and SIXJ, and the recoupling matrix is investigated by RE- 
COUPO (see section SAI_RECLS). 
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Table 1 : Scheme of the expressions for matrix elements of two-particle scalar operator (like Coulomb inter- 
action) . 
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(34) 


fO A \ 

(34) 


/o rr \ 

(35) 




(53) 


(26) 


(42) 


a/377 (52) 


(7) 


/O A \ 

(34) 


fO A \ 

(34) 


/o tr \ 

(35) 




(53) 


(26) 


(42) 


/3a77 (52) 


(7) 


/O A \ 

(34) 


fO A \ 

(34) 


(35) 




(53) 


(26) 


/ A OA 

(42) 


NONRELAT51 


















a/3~f5 (52) 


(8) 


/O A \ 

(34) 


/O A \ 

(34) 


(34) 


/O A \ 

(34) 


(53) 


(33) 


(43) 


/3a 7 (5 (52) 


(8) 


/O A \ 

(34) 


/a A \ 

(34) 


/O A \ 

(34) 


/O A \ 

(34) 


(53) 


(33) 


(43) 


a/3<57 (52) 


(8) 


/O A \ 

(34) 


/Q A \ 

(34) 


(34) 


fO A \ 

(34) 


(53) 


(33) 


(43) 


(3aSj (52) 


(8) 


/O A \ 

(34) 


/Q A \ 

(34) 


(34) 


/O A \ 

(34) 


(53) 


(33) 


(43) 


7<5a/3 (52) 


(8) 


/O A \ 

(34) 


(O A \ 

(34) 


(34) 


/O A \ 

(34) 


(53) 


(33) 


/ A 0\ 

(43) 


jSf3a (52) 


(8) 


/O A \ 

(34) 


(34) 


/O A \ 

(34) 


(34) 


(53) 


(33) 


(43) 


J 7 a/3 (52) 


(8) 


/O A \ 

(34) 


(34) 


(34) 


{ o A\ 

(34) 


(53) 


(33) 


(43) 


<5 7 /3a (52) 


(8) 


/O A \ 

(34) 


/Q A \ 

(34) 


/ o A\ 

(34) 


(34) 


(53) 


(33) 


/ A 0\ 

(43) 


NONRELAT52 


















a-y(35 (50) 


(8) 


/O A \ 

(34) 


/O /I A 

(34) 


/O /I \ 

(34) 


/O A \ 

(34) 


/cr 1 \ 

(51) 


(33) 


/ A 1 A 

(43) 


a 7 <5/3 (54) 


(8) 


fO A \ 

(34) 


(O A \ 

(34) 


/O A \ 

(34) 


/o /I \ 

(34) 


(55) 


(33) 


/ A O \ 

(43) 


7 a<5/3 (50) 


(8) 


/O A \ 

(34) 


/Q A \ 

(34) 


(34) 


fO A \ 

(34) 


/rr 1 \ 

(51) 


(33) 


(43) 


7 a/3(5 (54) 


(8) 


fO A \ 

(34) 


/O /I \ 

(34) 


(34) 


/O /I \ 

(34) 


/ rr rr \ 

(55) 


(33) 


/ A 1 A 

(43) 


(35aj (50) 


(8) 


/O A \ 

(34) 


{ o A\ 

(34) 


(34) 


fO A \ 

(34) 


/rr 1 \ 

(51) 


(33) 


(43) 


<5/3 7 a (50) 


(8) 


/O A \ 

(34) 


(34) 


/O A \ 

(34) 


(O A \ 

(34) 


/tr 1 \ 

(51) 


(33) 


(43) 


/3<57« (54) 


(8) 


/O A \ 

(34) 


/O A \ 

(34) 


(34) 


(34) 


/ rr rr \ 

(55) 


(33) 


(43) 


6f3aj (54) 


(8) 


(34) 


(34) 


(34) 


(34) 


(00 ) 


(33) 


(43) 


NONRELAT53 


















a<5/3 7 (50) 


(8) 


(34) 


(34) 


(34) 


(34) 


(51) 


(33) 


(43) 


<5a7/3 (50) 


(8) 


(34) 


(34) 


(34) 


(34) 


(51) 


(33) 


(43) 


a<5 7 /3 (54) 


(8) 


(34) 


(34) 


(34) 


(34) 


(55) 


(33) 


(43) 


SaPj (54) 


(8) 


(34) 


(34) 


(34) 


(34) 


(55) 


(33) 


(43) 


/37a<5 (50) 


(8) 


(34) 


(34) 


(34) 


(34) 


(51) 


(33) 


(43) 


7/3<5a (50) 


(8) 


(34) 


(34) 


(34) 


(34) 


(51) 


(33) 


(43) 


/3 7 <5a (54) 


(8) 


(34) 


(34) 


(34) 


(34) 


(55) 


(33) 


(43) 


7/3a6 (54) 


(8) 


(34) 


(34) 


(34) 


(34) 


(55) 


(33) 


(43) 
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For the distributions a[3otf3 and fiafia the subroutine NONRELAT1 uses (6) of [6], keeping in mind 
that 6 (n a \ a ,np\p 7 ~,) is expressed as (51) of [6] and tcnsorial parts B ( - Kl2ai ^ (n a X a , S), cK 2 ^ 2 ) (np\p, S) 
are equal to (35) from [6]. The coefficients 9 {n a X ai npXp, S) are investigated by COULOMBLS, the 

coefficients £( Kl2<712 ) (n a A Q , 5) and c( K ' 12<7 ' 12 ) (n^A^, S) are found by W1W2LS from the SAI_SQ library [5], 
and the recoupling matrix is calculated by RECOUP2. 

For the distributions a(3f3a and (3aaf3 the subroutine NONRELAT1 uses (6) of [6], keeping in mind 

that 6 (n a X a , npXp, S) is expressed as (55) of [6] and tcnsorial parts B (Kl2<Tl2) (n a X a , H), c( K ' 12 "' 12 ) (npXp, S) 
are equal to (35) from [6]. The coefficients (n a X a , npXp, S) are investigated by COULOMBLS and 

SIXJ, the coefficients B^ 2 ^ 2 ) (n Q A Q , S) and C^ 2 ^ 2 ) (n^, H) are found by W1W2LS from the SAI_SQ 
library [5], and the recoupling matrix is calculated by RECOUP2. 

COULOMBLS Investigates the two-electron submatrix elements of electrostatic interaction 

(niXinjXjllg^l^lln'iX'in'jX'^ , (34) 

according to the formula (9) of [7] . The values of these matrix elements are needed because of (see [6] ) 

(S) ~ (niXiHjXj \\g\\ n v X v nyXy) . (35) 

The value of the output parameter AA of this subroutine is: 

AA = 2[fc] 1 / 2 (/ i ||CW||^)(/ j ||C( fe )|^). (36) 

3.3 SAI_DUDU 

This library is meant for the calculation of matrix elements of any one- or two-particle operator. It contains 
44 subroutines. Most of the subroutines from this module use common blocks CONSTS and MEDEFN from 
package ATSP_MCHF [2, 3]. Similar to the SAIJJORE library it uses the methodology described in paper [6]. 
Therefore the arrangement of library SAIJXJDU is analogous to that of library SAIJJORE. Therefore we will 
not go into computational details of spin-angular parts of one- or two-particle operator matrix elements, 
but instead will reserve all the attention to a demonstration of connection of the Breit-Pauli operators to 
the general algorithm of one- or two-particle operator calculation. 

3.3.1 Spin own orbit interaction 

The subroutine SPINOR investigates the submatrix elements of the spin-own-orbit interaction operator 

(radl/i^lK-A;), (37) 

according to the formula (5) of [11]. 

3.3.2 Spin spin interaction 

The subroutine SSC investigates the spin-spin operator, which has the tensorial form (11) and finds the 
submatrix element of this operator between functions with any number of open shells. This subroutine is 
used for all distributions except aaaa, afiafi, fiafia, afifia, f3aaf3 and aa.fi (3. In the latter cases, instead 
of subroutine SSC the subroutines SS1111, SS1212, SS1212, SS1221, SS1221 and SS1122 are used. 

The subroutine SSA investigates the submatrix elements of the spin-spin interaction operator 

(n 4 A l n,A,||ffir i ^ 12 ^ 2 )||<A^;A;) (38) 

and 

(jiiXirijXj 1 1^-^+12,112) | K x >. n >. A / ) (39) 

respectively. This subroutine is used for all distributions except aaaa, afiafi, and fiafia. In the latter cases, 
instead of subroutine SSA the subroutine SSI is used. 
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3.3.3 Spin other orbit interaction 

The subroutine SOOC investigates the spin-other-orbit operator (16). This subroutine is used for all 
distributions except aaaa, a(3a(3, /3af3a, af3(3a, (5aaf3 and aa(3(i. For distributions aaaa, a(3a(3 and 
(3a(3a, according to (34), (39) and (40) of [12], the tensorial form 

ttsoo _ / rr(fc-lfel,101) , rj-(fc-lfcl,011) , rr(fc+lfcl,101) , rr(fe+lfel,011) \ 
n \2 — 2_^i\ n soo "i" n soo + a soo "T n soo J 

k 

is valid. Therefore instead of subroutine SSC the subroutines SOOllll, S001212, S001212 respectively 
are used. The distribution af3(3a, (3aaf3 and aa(3(3 are also calculated according to separate expressions. 
These are treated by the subroutines S001221P, S001221P and SOOH22P. 

The subroutines SOOA and SOOB investigate the submatrix elements of the spin-other-orbit interac- 
tion operator 

(nAnAI|tfL fe - lfel ' 101) lKA^.) , (nAnAI|tfL fe - lfel ' 011) lKA^.) , 

(^A^A.II^+^'^H^AXa;), (n i A i n j A j || J ff( s fe + lfel ' 011 )||^A^;A;.) 

according to the formulas (26), (27) and (28) of [12] respectively. These subroutines are used for all distri- 
butions except aaaa. In the latter case instead of subroutines SOOA and SOOB the subroutine SOOl is 
used. 



4 The Library usage for other programs 

There arc different versions of ATSP_MCHF [2, 3, 17], which support non-relativistic atomic structure calcula- 
tions. Each of them has some specific disparity. But it is possible to implement the library SAI in all of them 
and in other programs as well. The one example of such connection between ATSP_MCHF and SAI is shown 
for the program Breit from [3] in Table 2. The modify package ATSP_MCHF [3] which is fit to the module 
SAI is coming with the distribution of the library SAI. 

For writing the similar interface for other programmes like [18] we need to fill in the Common block 
MEDEFN from ATSP_MCHF, to write similar program as is presented in the Table 2 and to make some other 
subroutines similar as is presented in this distribution (the files breit. f, savenon.f and setupgg.f in the 
directory /dudu/breit). 

Table 2. THE MODIFY PROGRAM BREIT FOR CONNECTION 
BETWEEN ATSPJ4CHF AND SAI. 

PROGRAM BREITF 

IMPLICIT DOUBLE PRECISION CA-H,0-Z) 
PARAMETER (NWD=60 , NCD=1000 , NCD2=2*NCD) 
CHARACTER ANS*2, NAME(2)*24 

C0MM0N/DEBUG/IBUG1 , IBUG2 , IBUG3 , NBUG6 , NBUG7 , I FULL 

C0MM0N/C0NSTS/ZER0 , TENTH , HALF , ONE , TWO , THREE , FOUR, SEVEN , ELEVEN , EPS 
COMMON/DIAGNL/IDIAG , JA , JB 

COMMON /F0UT/N0V (2) , I0VLAP (10 , 2) , NF , NG , NR, NL , NZ , NN , NV , NS , I FLAG , NI J 
C0MM0N/IMAGNT/ CONST , C0NS00 , CONSS , ISPORB , IS00RB , ISPSPN , 
: IREL , I STRICT , IZOUT , IELST , ITENPR 

COMMON/INFORM/ IREAD , IWRITE , I0UT , ISC (8) 

C0MM0N/MEDEFN/IHSH,NJ(16) ,LJ(16) ,N0SH1(16) ,N0SH2(16) , J1QN1 (31 ,3) , 
: JiqN2(31,3) ,IJFUL(16) 

C0MM0N/PHASES/SIGNFA(NCD2) , ICSTAS 

COMMON/STATES/NCFG , NOCCSH (NCD2) , N0C0RB ( 5 , NCD2) , NELCSH ( 5 , NCD2) , 
: J1QNRD( 9,NCD2) ,MAX0RB,NJC0MP(NWD) ,LJC0MP(NWD) , IAJCMP (NWD) 

COMMON /OPERAT/ IC0L0M, IS0T0P , I0RB0RB 
COMMON /BREIT/ ISPORBG , IS00RBG , ISPSPNG 
COMMON/STEGG/IX , IGGG , IRHO , ISIG , IRHOP , ISIGP 
DIMENSION NFLGC20) ,IRFST(NCD2) ,NC0UNT(8) 
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EQUIVALENCE (NCOUNT(l) ,NF) 
LOGICAL INCL 

105 FORMAT (49H ISPORB=0 AND IS00RB=1 CAUSES THE PROGRAM TO FAIL, 
: 34H BECAUSE THE BLUME WATSON FORMULAE , / 
: 50H CANNOT BE USED FOR CLOSED SUBSHELLS . TO OVERCOME, 
: 34H THIS, THE CODE HAS SET ISPORB = 1//) 
11 FORMATC////' THE TYPE OF CALCULATION IS DEFINED BY ', 
: 'THE FOLLOWING PARAMETERS - '/ 

: 5X,22H BREIT-PAULI OPERATORS , 13X, 8HIREL =,12/ 

: 5X,27H PHASE CONVENTION PARAMETER, 8X.8HICSTAS =,12/) 
13 FORMAT (40H RELATIVISTIC OPERATORS TO BE INCLUDED -/5X,13H SPIN-ORB 

:IT (,I1,22H), SPIN-OTHER-ORBIT (,I1,15H), SPIN-SPIN (,I1,1H)/) 
24 FORMAT (36H0INITIAL DEBUG: IN 1 -ELECTRON PART =,I2,2H ,,5X, 

: 20H IN 2-ELECTRON PART=,I2/16X,23HIN RECOUPLING PACKAGE =,12/) 

42 FORMATC//' MATRIX ELEMENTS CONSTRUCTED USING ', 

: 'THE SPHERICAL HARMONIC PHASE CONVENTION OF'/) 

43 FORMAT (16X,47HC0ND0N AND SHORTLEY, THEORY OF ATOMIC STRUCTURE/ 16X , 
:47H ///) 

44 FORMAT (25X.42HFAN0 AND RACAH , IRREDUCIBLE TENSORIAL SETS/25X.42H— 
. ///) 

50 FORMAT (/20X, '=================================</ 

: 20X, 'BREIT-PAULI 2001'/ 

: 20X , ' ================================= • /) 

78 FORMAT (19H DEBUG PARAMETERS -/5X,16H NBUG6 (TENSOR) =,I2/5X, 
: ' NBUG7 (RELATIVISTIC OPERATORS - SO.SOO.SS) =>,I2//) 

* 

* SET INPUT AND OUTPUT CHANNELS 

ISOTOP=0 
IC0L0M=O 
IORBORB=0 
IREAD=4 
IWRITE=6 
IOUT = 7 
DO 2 I = 1,8 

ISC(I) =10+1 

NCOUNT(I) = 

OPEN(UNIT=ISC(I) ,STATUS=' SCRATCH' , FORM= ' UNFORMATTED ' ) 
2 CONTINUE 

NAME(l) = 'cfg.inp' 
NAME (2) = 'int. 1st' 

OPEN (UNIT=IOUT , FILE=NAME(2) , STATUS= 'UNKNOWN ' ) 

NIJ = 

NHDEL=10 

MXIHSH=(10) 
10 WRITE(IWRITE,50) 

WRITECO, ' (A/A/A/A) ') ' Indicate the type of calculation ', 
: ' => non-relativistic Hamiltonian only; ' , 
: ' 1 => one or more relativistic operators only;', 
: ' 2 => non-relativistic operators and selected relativistic: ' 

READ(5,*) IREL 

IF (IREL.NE.l) IC0L0M=1 

WRITE (0 , ' (A) ' ) ' Is full print-out requested? (Y/N) ' 
READ (5, ' (A2) ') ANS 
I FULL = 

IF (ANS . EQ. 'Y' .OR. ANS . EQ. 'y') I FULL = 1 
WRITE(0, ' (/A) ') 

: ' Phases:- Condon and Shortley or Fano and Racah ? (CS/FR) ' 
READ (5, ' (A2) ') ANS 

IF ( ANS . EQ. 'CS' .OR. ANS .Eq. 'cs') THEN 

ICSTAS = 1 
ELSE 

ICSTAS = 
END IF 
ISPORB = 
ISOORB = 
ISPSPN = 
IORBORB=0 

IF (IREL .NE. 0) THEN 
ISPORB = 1 
ISOORB = 1 
ISPSPN = 1 

WRITE(0 , ' (A) ' ) ' All relativistic operators ? (Y/N) ' 
READ (5, ' (A2) ') ANS 

IF ( ANS .Eq. 'N ' .OR. ANS .Eq. 'n ') THEN 
WRITECO, '(A)') ' Spin-orbit ? (Y/N) ' 
READ (5, ' (A2) ') ANS 

IF ( ANS .Eq. 'N ' .OR. ANS .EQ. 'n ') ISPORB = 
WRITE(0, ' (A) ') ' Spin-other-orbit ? (Y/N) ' 
READ (5, ' (A2) ') ANS 

IF ( ANS .EQ. 'N ' .OR. ANS .EQ. 'n ') ISOORB = 
WRITECO, '(A)') ' Spin-spin ? CY/N) ' 
READ (5, ' (A2) ') ANS 

IF ( ANS . Eq. 'N ' .OR. ANS .EQ. 'n ') ISPSPN = 
END IF 
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IF (ISPORB . EQ . . AMD . ISOORB . NE . 0) THEN 
WRITE (IWRITE.105) 

END IF 
END IF 
IBUG1 = 
IBUG2 = 
IBUG3 = 
NBUG6 = 
NBUG7 = 

WRITE (IWRITE, 11) IREL,ICSTAS 
ISPORBG = ISPORB 
ISOORBG = ISOORB 
ISPSPNG = ISPSPN 
CALL FACTRLC32) 

* 

* — READ IN THE SET OF CONFIGURATIONS 

CALL ACNFIGCNAME(l)) 
IDG = 
NZERO = NCFG 
NEW = NCFG 
DO 554 I = 1, NZERO 
554 IRFST(I) = 1 
ISTART = 1 

WRITECO , ' CA) ' ) ' All Interactions? (Y/N) : ' 
READ (5, ' (A2) ') ANS 

IF (ANS .NE. 'Y' .AND. ANS .NE. >y») THEN 
WRITECO , ' CA, 13 , A/A) ' ) ' Of the ',NCFG, 
: ' configurations, how many at the end are new? ', 

: ' How many configurations define the zero-order set?' 

READ C5,*) NEW, NZERO 
IF CNEW .EQ. 0) NEW = NCFG 

ISTART = NCFG - NEW + 1 

IF (NZERO .eq. 0) NZERO = NCFG 

IF (IREL .NE. .AND. NZERO .NE. NCFG) THEN 

WRITECO, '(A)') ' Rel interaction with all the zero block? 
READ C5, ' CA2) ') ANS 

IF CANS .NE. 'Y' .AND. ANS .NE. 'y') THEN 
DO 553 I = 1,(20) 
553 NFLGCI) = 

WRITECO,*) ' Define your reference set : F0RMC20I3)' 

READ C5,'C20I3)>) CNFLGCD ,1=1, C20)) 

DO 551 I = 1, CISTART-1) 

IRFSTCI) = 

DO 552 J = 1, C20) 

IF CNFLGCJ) .EQ. I) THEN 

IRFSTCI) = 1 

GO TO 552 
ENDIF 

552 CONTINUE 
551 CONTINUE 
ENDIF 

WRITECO,*) ' Diagonal rel corrections ? CY/N) : ' 
READ C5, ' CA2) ') ANS 

IF CANS.EQ. 'Y' .OR.ANS.EQ. 'y') IDG = 1 
ENDIF 

ISTRICT = 1 

WRITECO,*) ' Restricted Two-body interactions? CY/N) ; ' 
READ (5, ' CA2) ' ) ANS 

IF (ANS .NE. 'Y' .and. ANS .NE. 'y' ) ISTRICT = 
END IF 

* 

* ... Start the calculation 

DO 6 JA = ISTART, NCFG 
ICOUNT = 

* 

* ... Set up defining quantum numbers for each matrix element. 

CALL SETUP (JA,JB) 

IFCIBUG1.NE.0 .OR. IBUG2.NE.0) CALL VIJOUTCJA, JB) 
IF(IHSH.GT.MXIHSH) STOP 

* 

* ... TEST ON POSSIBLE RECOUPLING ORTHOGONALITY. 

CALL ORTHOG(LET, INCL) 
IF (LET .NE. 0) THEN 
CALL SETUPGG 
IF (IX .LE. 4) THEN 

IF CIFULL .NE. 0) WRITECIWRITE,77) 
77 F0RMATC///30X, 'MULTIPLYING FACTOR' , 1 IX, 'TYPE OF INTEGRAL') 

CALL NONBP 

IF CIFLAG .NE. 0) NIJ = NIJ + 1 
ENDIF 
ENDIF 
7 CONTINUE 
6 CONTINUE 

NTOTAL = CC2*NCFG - NEW +l)*NEW)/2 
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WRITE (0,*) NTOTAL , ' matrix elements' 

WRITECO,*) NIJ, ' non-zero matrix elements' 

WRITECO,*) 100*NIJ/REAL (NTOTAL) , ' 7. dense' 

WRITECO,*) 'NF=',nf,' KG=',ng,' NR=',nr,' NL=' ,nl 

WRITECO,*) 'NZ=',nz,' NN=',nn,' NV=',nv,' NS=' ,ns 

WRITECO,*) 'Total number of terms =', NF+NG+NR+NL+NZ+NN+NV+NS 

CALL OUTLSJ 

WRITE (IWRITE, 42) 

IF(ICSTAS.EQ.O) THEN 

WRITE(IWRITE,44) 
ELSE 

WRITE(IWRITE,43) 
END IF 
STOP 
END 



5 Installation of the Library 

The library SAI will be distributed as compressed tar archive file. On a UNIX (or compatible) workstation, 
the commands uncompress targ_lib_SAI and tar xvf targ_lib_SAI reconstruct the file structure. The 
directory LIBRARY _SAI contains the source code (subdirectories dudu, nore, reels, sqlsf 1 and sqlsf 1) as 
well as makefile make_lib_SAI. It also includes a modify package ATSP_MCHF [3] which is fited to the module 
SAI. 

The command make_lib_SAI creates the libraries libdudu . a, libnore . a, librecls . a, libsqlsf 1 . a and 
libsqlsf 2 . a in directory bin. 

6 Examples 

The way to determine the ionization potential for Ce is demonstrated in this section. For this, one needs 
to calculate the ground state energy of Ce atom ( 1 G), the ground state of a Ce+ ion ( 2 G), and find the 
difference. It has been shown in [19, 20, 21] that the non-relativistic Hartree-Fock ionization potential 
of singly ionized lanthanides follows the observed trends. Therefore we have chosen the non-relativistic 
approximation, and find the energy values by solving the multi-configuration Hartree-Fock equations. For 
brevity, we present only the Ce + test cases. The energy for Ce is found in the same way. 

In Table 3 we present an example of generating a set of configurations where the ground state energy 
of th Ce + ion is to be calculated by the multiconfiguration Hartree-Fock method. The configurations are 
obtained by taking the singles and doubles excitations from {4/, 5d, 6s} to the active set {5/, 5g} (since 
neither 5d nor 6s are filled shells, better accuracy would be obtained if these orbitals were also included in 
the active set) [21]. In this case the program GENCL [22] generates a cfg.out file with 40 configuration states. 

Table 3. TEST RUN OUTPUT 



»genclf 



You are under the program GENCL 

which GENerates a Configuration List 
Type H CHelp) to obtain information about the input format 
Type < RETURN > if you already know 



Header ? 

> Ce+ ground state 

Closed Shells ? 

> Is 2s 2p 3s 3p 3d 4s 4p 4d 5s 5p 

Reference Set ? 

> 4f Cl)5d(l)6sCD 

2 ? 

> 

Active Set ? 

> 

Replacements ? 

> sd 

Virtual Set ? 

> 5f,5g 

From which shell ? 

> 1 
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> 3 



> 2G 



To which shell ? 



Final Terms ? 



INPUT D 
Ce+ ground sate 
Is 2s 2p 3s 
4f (l)Bd(l)6s(l) 



3p 3d 4s 4p 4d 5s 5p 



Closed Shells 
Configuration 
Their couplings 



Closed Shells 
Reference set 
Virtual Set 
S-Replacement 
Configuration 
Their couplings 



Closed Shells 
Reference set 
Virtual Set 
S-Replacement 
Configuration 
Their couplings 



5s 5p 



3d 4s 4p 4d 5s 5p 



*************** INPUT DATA ********** 

Header 
Closed shell 
Reference Set 
Active Set 
Replacements : sd 
Virtual Set : 5f,5g 

From which shell : 1 
To which shell : 3 
Final Terms : 2G 
GENERATE ALL COUPLINGS FOR EACH MEMBER OF THE REFERENCE SET 

Ce+ ground state 

Is 2s 2p 3s 3p 3d 4s 4p 4d 
4f( 1) 5d( 1) 6s( 1) 

2F1 2D1 2S1 1G0 2G0 

2F1 2D1 2S1 3G0 2G0 

FOR VIRTUAL SET GENERATE CONFIGURATION STATES FOR S-REPLACEMENT 

Ce+ ground state 

Is 2s 2p 3s 3p 
4f (l)5d(l)6s(l) 
5f ,5g 
4f = 5f 

5d( 1) 6s( 1) 5f( 1) 

2D1 2S1 2F1 1D0 2G0 

2D1 2S1 2F1 3D0 2G0 

Output omitted for brevity . . . 

FOR VIRTUAL SET, GENERATE CONFIGURATION STATES FOR D-REPLACEMENT 

Ce+ ground state 

Is 2s 2p 3s 3p 
4f (l)5d(l)6s(l) 
5f ,5g 

4f .5d =5f .5g 
6s( 1) 5f( 1) 5g( 1) 

2S1 2F1 2G1 1F0 2G0 

2S1 2F1 2G1 3F0 2G0 

Output omitted for brevity . . . 
OK! 

List of configurations and their couplings 
is in the file cfg.inp 



3d 4s 4p 4d 5s 5p 



The generating of the angles file int.lst is demonstrated in Table 4 with the modify package ATSP_MCHF [3] 
which is fited to the module SAI and presented in the distribution. It will be needed for calculating the 
energy of 2 G state of Ce+ with the MCHF program [23]. The program input file cfg.inp is taken with the 
help of renaming cfg.out to cfg.inp. 



Table 4. TEST RUN OUTPUT 

>>breitf ========> Obtain Expressions for 

Breit-Pauli Interactions 



BREIT-PAULI 



Indicate the type of calculation 

=> non-relativistic Hamiltonian only; 

1 => one or more relativistic operators only; 

2 => non-relativistic operators and selected relativistic: 
>0 

Is full print-out requested? CY/N) 
>N 

Phases:- Condon and Shortley of Fano and Racah ? (CS/RF) 
>CS 

All Interactions? (Y/N) 
>Y 

THE TYPE OF CALCULATION IS DEFINED BY THE FOLLOWING PARAMETERS - 
BREIT-PAULI OPERATORS IREL = 

PHASE CONVENTION PARAMETER ICSTAS = 1 



THE CONFIGURATION SET 
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STATE (WITH 40 CONFIGURATIONS) : 



THERE ARE 5 ORBITALS AS FOLLOWS: 

4f 5d 6s 5f 5g 
THERE ARE 11 CLOSED SUBSHELLS COMMON TO ALL CONFIGURATIONS AS FOLLOWS: 

Is 2s 2p 3s 3p 3d 4s 4p 4d 5s 5p 

CONFIGURATION 1 ( OCCUPIED ORBITALS= 3 ): 4f ( 1) 5d( 1) 6s( 1) 
COUPLING SCHEME: 2F1 2D1 2S1 

1G0 2G0 

. . . Output omitted for brevity . . . 

CONFIGURATION 40 ( OCCUPIED ORBITALS= 2 ): 4f( 1) 5g( 2) 
COUPLING SCHEME: 2F1 3H2 

2G0 

MATRIX ELEMENTS CONSTRUCTED USING THE SPHERICAL HARMONIC PHASE CONVENTION OF 
CONDON AND SHORTLEY, THEORY OF ATOMIC STRUCTURE 



»cat int. 1st ===> Display the int. 1st file produced 

Ce+ ground sate 



F 


C 


4f , 


6d) 


12 


R 


( 


4f 


5d, 


5f 


5d) 


14 


R 5( 


5d 5g, 5f 


5f ) 


2008 


F 


C 


4f , 


6s) 


16 


R 


C 


4f 


6s , 


5f 


6s) 


20 


R 5( 


5f 5f , 5g 


5g) 


2014 


F 


( 


4f , 


5f) 


22 


R 


OC 


4f 


5g> 


5f 


5 g) 


32 


R 6( 


4f 5d, 5f 


5g) 


2018 


F 


0( 


4f , 


5g) 


40 


R 


1 c 


4f 


5d, 


5d 


5f ) 


46 


R 6( 


4f 5g , 4f 


5d) 


2082 


F 


0( 


5d, 


6s) 


44 


R 


1 c 


4f 


5d, 


5g 


5f ) 


50 


R 6( 


4f 5g, 5f 


5d) 


2146 


F 


0( 


5d, 


5f) 


56 


R 


1( 


4f 


5f , 


5d 


Bg) 


110 


R 6( 


4f 5g, 5f 


Bg) 


2178 


F 


0( 


5d, 


Bg) 


76 


R 


1( 


4f 


5f , 


5g 


5d) 


124 


R 6( 


5d 5f , 5g 


5f ) 


2182 


F 


0( 


6s, 


5f) 


80 


R 


1 ( 


4f 


5g> 


5d 


4f ) 


188 


R 6( 


Sd 5g, 5g 


Bg) 


2232 


F 


0( 


6s, 


5g) 


84 


R 


1 c 


4f 


Sg> 


5d 


5f ) 


250 


R 7( 


4f 5g, 5g 


5f ) 


2336 


F 


0( 


5f , 


5f) 


90 


R 


1 ( 


4f 


5g. 


5g 


Bf ) 


354 


R 7( 


5f 5f , 5g 


5g) 


2342 


F 


0( 


5f , 


5g) 


102 


R 


1 c 


5d 


5f , 


5f 


Bg' 


356 


* 








F 


0( 


Bg. 


5g) 


108 


R 


1 ( 


5d 


5g> 


5f 


Bf ) 


406 


-1 


OOOOOOOOR 


26 


14 


F 


2( 


4f , 


5d) 


120 


R 


1 c 


5f 


5f , 


5g 


Bg) 


412 


1 


OOOOOOOOR 


27 


15 


F 


2( 


4f , 


5f) 


130 


R 


2 C 


4f 


5d, 


4f 


6s) 


414 


1 


OOOOOOOOR 


23 


11 


F 


2( 


4f , 


5g) 


166 


R 


2 C 


4f 


5d, 


5f 


5d) 


428 





50000000R 


3 


1 


F 


2( 


5d, 


5f) 


178 


R 


2 C 


4f 


5d, 


5f 


6s) 


430 





86602540R 


4 


1 


F 


2( 


5d, 


5g) 


226 


R 


2 C 


4f 


5d, 


5f 


Bg) 


434 


-1 


OOOOOOOOR 


22 


10 


F 


2( 


5f , 


5f) 


232 


R 


2 ( 


4f 


6s , 


5f 


5d) 


436 


1 


OOOOOOOOR 


21 


9 


F 


2( 


5f , 


5g) 


258 


R 


2 C 


4f 


5 g> 


4f 


5d) 


500 






19 


■j 


F 


2( 


5g, 


5g) 


264 


R 


2( 


4f 


5g, 


5f 


5d) 


564 


-1 


OOOOOOOOR 


20 


8 


F 


4( 


4f , 


5d) 


276 


R 


2( 


4f 


Bg, 


5f 


Bg) 


604 


-1 


OOOOOOOOR 


28 


16 


F 


4( 


4f , 


5f) 


287 


R 


2( 


5d 


5f , 


6s 


Bf) 


606 


1 


OOOOOOOOR 


25 


13 


F 


4( 


4f , 


5g) 


324 


R 


2( 


5d 


5f , 


5g 


5f) 


610 


-0 


86602540R 


3 


2 


F 


4( 


5d, 


5f) 


336 


R 


2( 


5d 


5g, 


6s 


5d) 


650 





50000000R 


4 


2 


F 


4( 


5d, 


5g) 


388 


R 


2( 


5d 


5g, 


6s 


Bg) 


670 


-1 


OOOOOOOOR 


24 


12 


F 


4( 


5f , 


5f) 


394 


R 


2( 


5d 


5g, 


5g 


Bg) 


720 


-0 


86602540R 


3 


2 


F 


4( 


5f , 


5g) 


420 


R 


3( 


4f 


5d, 


5d 


5f) 


734 





50000000R 


3 


1 


F 


4( 


5g, 


Bg) 


426 


R 


3( 


4f 


5d, 


6s 


4f) 


736 





86602540R 


4 


1 


F 


6( 


4f , 


5f) 


435 


R 


3( 


4f 


5d, 


6s 


5f) 


738 


-1 


OOOOOOOOR 


18 


6 


F 


6( 


4f , 


5g) 


464 


R 


3( 


4f 


Sd, 


5g 


5f) 


742 





50000000R 


4 


2 


F 


6( 


5f , 


5f) 


470 


R 


3( 


4f 


6s, 


5d 


5f) 


744 


1 


OOOOOOOOR 


17 


5 


F 


6( 


5f , 


5g) 


490 


R 


3( 


4f 


6s, 


6s 


5f) 


748 


-1 


OOOOOOOOR 


18 


6 


F 


6( 


5g, 


5g) 


496 


R 


3( 


4f 


6s, 


5g 


5f) 


768 


1 


OOOOOOOOR 


25 


13 


F 


8( 


Bg. 


5g) 


502 


R 


3( 


4f 


5f , 


5d 


6s) 


780 


1 


OOOOOOOOR 


17 


5 


* 










R 


3( 


4f 


Sf , 


5d 


Bg) 


836 


1 


OOOOOOOOR 


19 


7 






OOOOOOOOF 


1 1 


R 


3( 


4f 


5f , 


6s 


5d) 


848 


-1 


OOOOOOOOR 


20 


8 






OOOOOOOOF 


13 13 


R 


3( 


4f 


5f , 


6s 


Bg) 


860 


-1 


OOOOOOOOR 


24 


12 






OOOOOOOOF 


14 14 


R 


3( 


4f 


5f , 


5g 


5d) 


884 


1 


OOOOOOOOR 


21 


9 






OOOOOOOOF 


7 7 


R 


3( 


4f 


5f , 


5g 


6s) 


890 


-1 


OOOOOOOOR 


26 


14 






OOOOOOOOF 


15 15 


R 


3( 


4f 


Bg, 


5d 


4f) 


942 


-1 


OOOOOOOOR 


28 


16 


Ccont 


) 






R 


3( 


4f 


Bg, 


5d 


5f) 


992 


-1 


OOOOOOOOR 


22 


10 








* 




R 


3( 


4f 


Sg, 


6s 


4f) 


1024 


1 


OOOOOOOOR 


23 


11 


G 


0( 


4f , 


5f) 


21 


R 


3( 


4f 


Bg, 


6s 


5f) 


1046 


1 


OOOOOOOOR 


27 


15 


G 


1( 


4f , 


5d) 


33 


R 


3( 


4f 


Bg, 


5g 


5f) 


1150 





17142857R 


21 


9 


G 


1( 


4f , 


5g) 


111 


R 


3( 


5d 


6s, 


5f 


5f) 


1152 





02857143R 


19 


7 


G 


1( 


5d, 


5f) 


124 


R 


3( 


5d 


5f , 


5f 


6s) 


1154 





14285714R 


4 


2 


G 


1( 


5f , 


Bg) 


182 


R 


3( 


5d 


5f , 


5f 


Bg) 


1156 


-0 


14285714R 


3 


1 


G 


2( 


4f , 


5f) 


202 


R 


3( 


5d 


Bg, 


5f 


Bf) 


1206 





08571429R 


20 


8 


G 


2( 


5d, 


6s) 


207 


R 


3( 


6s 


5f , 


5f 


Bg) 


1216 





28571429R 


22 


10 


G 


2( 


5d, 


5g) 


317 


R 


3( 


6s 


Bg, 


5f 


5f) 


1218 





42857143R 


28 


16 


G 


3( 


4f , 


5d) 


329 


R 


3( 


5f 


5f , 


5g 


5g) 


1224 





28571429R 


27 


IB 


G 


3( 


4f , 


6s) 


334 


R 


4( 


4f 


5d, 


5f 


5d) 


1238 


-0 


24743583R 


4 


1 


G 


3( 


4f , 


5g) 


413 


R 


4( 


4f 


5d, 


5f 


Bg) 


1242 





02857143R 


24 


12 


G 


3( 


5d, 


5f) 


426 


R 


4( 


4f 


6s, 


5f 


Bg) 


1252 





42857143R 


23 


11 


G 


3( 


6s, 


5f) 


431 


R 


4( 


4f 


Bg, 


4f 


5d) 


1316 





17142857R 


26 


14 


G 


3( 


5f , 


5g) 


489 


R 


4( 


4f 


Bg, 


4f 


6s) 


1338 


-0 


24743583R 


3 


2 



1G 



G 4( 


41 , 


ol ) 


510 


R 




4f 


5 g. 


5f 


ouj 


1402 


. 




25 


13 


G 4( 


5d, 


5g) 


606 


R 


4( 


4f 


5 g. 


5f 


6s) 


1434 


-0 . 


14982984R 


17 


1 


G 4( 


6s , 


5g) 


612 


R 


4( 


4f 


5 g> 


5f 


5g) 


1476 


. 


Jb9ol2o9R 


17 


2 


Got, 


41 , 


5d) 


624 


R 


4( 


5d 


6s , 


5g 


5g) 


1478 


. 


14yo2yo4K 


18 


2 


G 5 C 


41 , 


5g) 


703 


R 


4C 


5d 


5f , 


5g 


5f ) 


1482 


. 


Jo9ol2cSyK 


18 


1 


G 5( 


5d, 


5f) 


716 


R 


4( 


5d 


5 g. 


5d 


6s) 


1512 


-0 , 


03937627R 


30 


28 


G 5( 


5f , 


5g) 


774 


R 


4C 


5d 


5 g. 


5g 


6s) 


1552 


-0 . 


. 08460778R 


32 


25 


G 6 C 


4f, 


5f) 


795 


R 


4 ( 


5d 


5 g. 


5g 


5g) 


1602 


. 


1 loooyUbR 


33 


24 


G 6 C 


5d, 


5g) 


905 


R 


4 C 


6s 


5f 


S g 


5f ) 


1622 


0. 


00067076R 


34 


28 


G 7( 


41 , 


5g) 


984 


R 


4( 


6s 


5g! 


5g 


Sg) 


1624 


-0. 


05155563R 


30 


27 


G 7( 


5f , 


5g) 1042 


R 


5( 


41 


5d, 


5d 


51) 


1638 


-0, 


04956348R 


33 


23 


* 








R 


5( 


41 


5d, 


5g 


5f) 


1642 


0. 


12698464R 


30 


26 


-0 


.27380952G 


34 34 


R 


5( 


41 


5f , 


5d 


5g) 


1702 


0. 


09770956R 


30 


25 


-0 


.50787450G 


34 33 


R 


5( 


41 


5f , 


5g 


5d) 


1726 


0. 


12978765R 


33 


22 





.55135887G 


34 32 


R 


5( 


41 


5g, 


5d 


41) 


1790 


-0. 


21161602R 


30 


24 


-0 


.09731237G 


34 31 


R 


5( 


41 


5g, 


5d 


5f) 


1852 


-0. 


. 12989160R 


33 


25 


-1 


.00946357G 


34 30 


R 


5( 


41 


Sg, 


5g 


5f) 


1956 


0. 


02273390R 


30 


23 


Ccont 


.) 






R 


5( 


5d 


51, 


5f 


5g) 


1958 


(cont . ) 







The energy of 2 G state of Ce + is found by the multi-configurational Hartree-Fock method. First, the 
wave functions are obtained by the HF method, using the HF96 program (see Gaigalas and Fischer [24]). 
These functions are used as starting ones in solving the multi-configurational Hartree-Fock equations. The 
input files for MCHF program [23] int.lst, and cfg.inp, are taken from the examples described above. 

7 Conclusion 

The library SAI supports large scale computations of open-shell atoms using multi-configuration Hartree- 
Fock or configuration interaction approaches and may be useful for developing codes for calculating the 
spin-angular parts of effective operators from many-body perturbation theory and orthogonal operators or 
for evaluating the rclativistic Hamiltonian in LS'-coupling as well as for various versions of semi-empirical 
methods. For example, it expands the use of possibilities of package MCHF_ASP [2, 3, 17] because the new 
program is faster and provides coefficients and matrix elements for shells (nl) with I — 0, 1, 2 and 3, and I 2 
for I > 3. It is important to underline that the use of the library described allows to perform calculations 
of energy spectra of practical any atom or ion starting with LS-coupling and accounting for all relativistic 
corrections of the order a 2 . 

Program is obtainable from the Institute of Theoretical Physics and Astronomy, A. Gostauto 12, Vilnius, 
2600, LITHUANIA. E-mail: gaigalas@itpa.lt. 
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